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Abstract 

We determine the energy spectrum and the corresponding eigenfunctions of a 2D Dirac 
oscillator in the presence of Aharonov-Bohm (AB) effect . It is shown that the energy 
spectrum depends on the spin of particle and the AB magnetic flux parameter. Finally, 
when the irregular solution occurs it is shown that the energy takes particular values. The 
nonrelativistic limit is also considered. 
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1 Introduction 



The Dirac oscillator was introduced for the first time by Moshinsky and Szczepaniak pQ. The 
basic idea is that the momentum p is replaced in the free Dirac equation by p — > p — iMuftr, 
where r is the position vector, M is the mass of particle and u) the frequency of the oscillator. 
The Dirac oscillator has attracted much attention and has various physical applications pSOH], 
particularly in semiconductor physics [5]. In the nonrelativistic limit, the Dirac oscillator 
becomes an harmonic oscillator with strong spin-orbit coupling term. 

Recently, the 2D electron systems has become an active research subject due to the advances 
in nanofabrication technology like quantum wells, quantum wires, quantum dots, quantum Hall 
effect and high superconductivity [HIE]. Currently, the 2D parabolic potential V = \Moj 2 r 2 is 
often used to describe confined 2D systems in nonrelativistic case. For relativistic case, such 
systems are described using the 2D Dirac oscillator [3]. 

On the other hand, the Aharonov-Bohm (AB) effect [7] leads to a number of remarkable 
interference phenomena in mesoscopic systems [8], and has been used also to study the 2D 
model in both superconductivity and particle theory [5]. In this article, in order to study 
a bound state version of the AB effect, let us consider the 2D Dirac oscillator of relativistic 
particle in the presence of the AB effect. The associate magnetic field is assumed to be confined 
to a filament of vanishingly small radius and perpendicular to the plane where the particles are 
confined to move. We study the dependence of the bound states energy and eigenf unctions on 
the AB magnetic flux parameter and the spin of particle. 

This paper is organized as follows. In Section 2, we solve the Dirac oscillator in polar 
coordinates with an additional delta-function interaction due to the AB effect. In Section 3, we 
analyze and discuss the energy spectrum and its degeneracy. The nonrelativistic limit is also 
considered. Along the paper we adopt the natural units, H = 1, c = 1 . 



2 Solution of the Dirac equation 

In the beginning, one writes the stationary Dirac equation for a particle of mass M in terms of 
two-component spinors ip 

Eil) = (/?7tt + (3M)ij, (1) 
where it = p — eA is the minimally coupling and A is the vector potential. 
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The 2D Dirac oscillator is defined by changing the momentum p in the Dirac equation 
[SI ID] as 

p -> p - iMcufJr. (2) 

The magnetic field H, associated to the AB effect, is assumed to be perpendicular to the plane 
and confined to a filament of vanishingly small radius so that the flux [HJ [TT] 



a = —e H(r)rdr, (3) 



o 



is finite and nonzero. Since, however, one would ultimately like to pass to the limit in which 
the field is confined, H must be taken to be of the form [21 CD] 

eH = ~ S (r-R), (4) 

the corresponding form of the potential A in the Coulomb gauge is [21 ED] 

f -2u- r>R, 
eA = ( r 5 
(0, r<R. 

Then, the stationary Dirac equation of a 2D Dirac oscillator in the presence of the AB effect is 

Eij) = [/3 7 (tt - iMu/3r) + (3M] ip. (6) 

Since we are using only two component spinors, the matrices (3 and /3ji are conveniently defined 
in terms of the Pauli spin matrices as 

Py l = a\ [3 l2 = sa\ (3 = a\ s = ±1. (7) 

Eq. © can be rewritten as 

P [7 (tt - iMuj(3r) + M — f3E] ip = 0. (8) 
By applying the following matrix operator 

[7 (tt - iMufir) -M-(3E}f3 

to Eq. (JH1), we obtain 

[7 (tt - iMupr) 7 (tt - iMu(3r) - M 2 + E 2 ]ifj = 0. (9) 
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The first term in the left hand side can be computed as 



7 (vr - iMuj(3r) 7 (tt - iMufir) = jit^n + 2Muf3 + 2Mus(n Ar)J - M 2 u 2 r 2 

= yiryx + 2Mu(3 - 2Mus [L z + a}(3- M 2 u 2 r 2 . 



Thus, Eq. © reads 



[ 7 7T7^ - M 2 u 2 r 2 + 2Mu [1 - s (L z + a)} (3 + E 2 - M 2 ] tj) = 0, 



(10) 



where 



Then, in the coordinates representation, Eq. (JTOl) reads 



77T77T = — 7T 2 + esa^H, L z 



ld_ 



-^r^+-(^-+ta) 2 -M 2 u J 2 r 2 +E 2 -M 2 +2Mu 
r or Or r z oa 



A 9 s 



(11) 



we note that the last equation contains a spin-orbit coupling term. 
Let us define the first component as 



irnO 



(12) 



m=— 00 



where the radial part satisfy the equation 



d 2 1 d 
dr 2 r dr 



_ , r — _ M Vr 2 - (m ^ 2 a) + £ 2 - M 2 + 2M^ [1 - s (m + a)} - ^6(r - R) 



Thus, it yields the following equations for the two regions 



d 2 Id , 9 m 2 , 9 

I A r |- k 

dr 2 r dr r 2 " 



+ -^-AV- (m + a)2 +fc 2 



r d 2 Id 

dr 2 r dr 



/ m (r) = 0, r<R, 
f m (r) = , r > i?, 



where 



fc 2 u< = E 2 - M 2 + 2Mu [1 - s (m + a)] , 
fc t ? n = E 2 — M 2 + 2Mw (1 - sm) , A 2 = MV. 



/m(r) = 0. 
(13) 

(14) 
(15) 

(16) 
(17) 
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The general solution of the radial equation (1151) is given in terms of the confluent hypergeometric 
function 



m + a + 1 



k 



out 



+ B m r- (m+a ^ 1 F 1 [ - 



1 — (m + a) — 



k 2 

^out 

2Mw 



2Mw j 
, 1 - (to + a) , Mur 2 



, m + a + 1, Mur' 



18) 



where A m , S m and C m are constants. 
This equation can be rewritten as 

/ m (r)=exp(-^r 2 )|yl m rl'"+'>l 



''■"''5 



|m + a| + 1 



A- 



out 



1 



1 — |m + at I 



out 

2Mu 



2Mcj j 
1 — Ito + a I , Muir 2 



m + a\ + 1, Mur' 



and consequently for the two regions we have 

/ m (r) m = C m rl m I exp(-^r 2 ) ^ Q 
for r < R, and 



Iml + 1 



',,2 



2Mw 



, |m| + 1, Mur' 



(19) 



(20) 



Jm{r)out = exp( — 



r) A m rl m+a l 



m + a| + 1 
1 — \m + a 



k 2 

^out 



m + a I + 1, Muir A 



h 2 

^out 



2Mu^ 
, 1 — |m + a| , Mwr 2 



(21) 



2Mw_ 
for r > R. 

The effect of the delta function is then taken into account by means of the continuity 
relations 

fm{R — £)in = fm{R + £)out , (22) 



dfmjr) 
dr 



R+e 



J R-e 



(23) 



In the lowest order in i? the confluent series and its first derivative behave like 



[xFx^c, Mur 2 )] r=R ^l, 



d x Fx (a, c, Mur 2 
dr 



2a 



MuR, 



r=R 



so Eq. (1221) gives 



(24) 
(25) 

(26) 



By applying Eq. (1221) and Eq. (I23p . one obtains for the unnormalized f m (r) for r > R (to 
lowest order in R) 



f m (r) = R lml <R 



\m+a\ 



\m\ + \m + ct\ + as + R 2 Mu 


( \m\ 


+ 1 ^ n ) (1 


m + a 


k 2 \ -, 

2Mu)> 




\m 


+ 1 1- 


m + a 





x r 



777 + a + 1 — 



2|m + a| + R 2 Mu 
h 2 



k 2 

2Mu 



k 2 

2Mu 



out 



+R\ 



x r 



2Mu 

(|m| + \m + a| + as) + R 2 Muj 



m + a\ + 1 1 — \m + a| 
, |m+ a | + 1, Mur 2 

kL 



(\ m \ + 1 -2Mz) (l-\m + a\- 



k 2 



2Mw 



\m\ + 1 



1 — m + a 



2 |m + a | + R 2 Mu 

A 



k 2 

2Mlu 



k 2 

"■out 

2Muj 



m+a+1 1— m + a 



1 — m + a — 



out 

2Mu 



1 — \m + a\ , Mur > exp( 



(27) 



By taking the limit i? — > 0, (i? 2 ~ 0), one clearly obtains 



771 + as 



2 777 + a 



„|m+a| 



x iFi 



x ^1- 



1 777 + a| + 1 — 



out 



777 + a 



2Mu 
h 2 

h out 

2Mu 



m 



J \2 2\m + a\J 



-|m+a| 



1 777 + a| , Mur' 



exp 



Mu 



(28) 



We note that the computed coefficients of the last equation are the same as in [9]. Thus, 
the solution of Eq. ffT3"]) in the limit R — > must always be the regular one unless Eqs. (9) and 
(10) of [9] are both satisfied. In this latter case only the irregular solution is allowed. 

Eqs. (11. a) and (ll.b) in [9] claim that the irregular solution occurs when 



(a) 777 = -N, N > 0, s = -1 or 

(b) m = -N-l, iV+l<0, s = +l, 



(29) 
(30) 
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then, the wave function ipx is given by 



+00 , 
Mr,0)= a mr lm+al iFA- 



\m + a\ + 1 



+ 6(-s)6(a)a„ N e- im r-t X F X 



1-e 



k 2 
k 2 



\m + a\ + 1, Mur e 2 r e 



2Mu 



+ d{s)6{-a)a. N . 1 e- l{N+1)e r^ 1 1 F 1 



, 1 -£,Mur z J e~~ 
,£,M^r 2 ) e 



2j g-^-'' J - '<«<-' 
2 \ MuiJl 



(31) 



where the prime on the summation indicates the omission of the two terms specified by 



Eqs. 



and 



a = N + £, < £ < 1 and iV is fixed integer. 



(32) 



The confluent series behaves asymptotically at large positive values of its argument as 



T(c) 

F(a,c,z) -> =r4e z z a - c 



r a) 



(33) 



so, the wave function tpi behaves as 



+00/ 



V"- | m+Q | M. r 2 -i(^+|m+a|+l 



+ ^(-s)0(«)a_ 7V e" iJV V^e— r r 
+ 6>(s)6>(— a)a_Ar_ie 



(34) 



where 6>(x) is the Heaviside step function and constant. 

The last expression is exponentially divergent. The divergence cannot be avoided except 
by putting the parameter [12J a = —n, with n = 0,1, 2, thus transforming the series into a 
polynomial of degree n (Laguerre polynomials). Hence 



Mw 2 

ur" e 2 e 



+ 9(-s)9(a)a„ N>n r^^(Mwr 2 )e-^ r V™ 
+ 0(s)0(-a)a_*-i,„ r^L^MwrV^e^"^, 

n\T(\m+a\+l) 
~ r(|m+o[+l+n) 0m ' 



where the constant a r 



(35) 
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When the conditions (I29I) - (I30I) are not satisfied, the energy spectrum for the regular solution 
is obtained from the condition 



1 

2 

thus, it follows from Eq. flHol) that 



\m + a\ + 1 



v out 

2Mu 



E = ±^M 2 + 2Mu [\m + a\ + s (m + a) + 2n}. 

The energy spectrum when the irregular solution occurs, for the first case (m 
s = —1), is deduced from 



(36) 



(37) 



-N, N > 0, 



1-^ 



k 2 
2Mu 



-n , 



(3? 



then 



For the second case (m 



E = ±^/M 2 +AMlo [ra-f]. 
-N- 1, iV + 1 < 0, s = +1) 



P 

2Mu 



-n, 



(39) 



(40) 



it follows that 



E = ±,JM 2 + AMcu [n + £ - 1]. 



(41) 



We remark that the energy levels given by Eqs. ( l39i) -( |4TI) depend only on the AB parameter £ 
defined by Eq. (JS2D. 

The lower component can be computed from Eq. 



Eifj = [/? 7 r 7T r . + (3-f\e ~ iMuf3frf3 + f3M] V, 



where 



Y = 7 1 cos 6 1 + 7 2 sin 9, 
*~f = — 7 1 sin + 7 2 cos 9. 
This leads to two first-order coupled differential equations 

Ei> 2 = e ise (ir r + isiro - iMujr) ^1 - M^ 2 



(42) 

(43) 
(44) 



(45) 



The second equation in the set of Eqs. (|4"5"]) leads to 

ip 2 = — — re"* [ — ~s --s- + Mur )tb 1 . 46 

r (E + M) \dr r i d9 r J Y K J 

Using the relation 

dL^Mur 2 ) = _ 2MurL , t i (aW) ; n > Xj (47) 
and substituting Eq. into Eq. (|4U|) one obtains 



• _M<£„.2 ^ +oo' 

-ze e 2 



+ M) 
+ 6(-s)9(a)a~ Nin 



m=—oo 



[fM-s(a + m)) r"- l L^(Mur 2 ) - 2Mwr" +1 L^\(Mur 2 



[-£ + s (TV - a)] r-^L-^Mur 2 ) - 2Mur^ +1 L^ 1 (M l 



ur 2 ) 



e -iN9 



+ 9(s)9{-a)a_ N _i :n 



2Mujr^L li n _ l (Mujr 



where /i— \m + a\ 



2\ 



[£ - 1 - s (N + 1) - aa] r^L^M^r 2 ) 
e -i(*+i)*\ (4f 



3 Analysis of energy spectrum 

Let us discuss the energy spectrum corresponding to the regular solution, i.e. 



E = ±^M 2 + 2Mu [\m + a\ + s (m + a) + 2n}. (49) 

From this relation, it is clear that the energy spectrum depends on the spin projection and the 
AB magnetic flux parameter. Two cases have to be distinguished: 
(a) For m + a > 0, Eq. (14*9"]) becomes 



E = ±^M 2 + 2Mu [(m + a) (s + 1) + 2n\. (50) 

For s = 1, it is to be noted that every state with (m, n ) has the same energy as that of a state 
with (m =p £,n± £) where £ is an integer. For s = —1, the energy depends only on the natural 
number n. 

(b) For m + a < 0, Eq. (|49p becomes 

E = ±^/M 2 + 2Mu [(m + a) (s - 1) + 2n\. (51) 
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For s = — 1, we observe that every state with (to, n ) has the same energy as that of a state 
with (to ± £, n ± £) where £ is an integer. For s — 1, the energy depends again only on the 
natural number n. 

On the other hand, we note that the energies Eqs. (|39l) - (l4Tl) corresponding to the two cases 
when the irregular solution occurs are not degenerated. 

3.1 Nonrelativistic limit 

Let us now examine the nonrelativistic limit of Eq. (11) by setting E = M + e, and considering 
e << M , so that E 2 » M 2 + 2M e, we obtain 

{^M + \ Mu?r2 ~ U[l ~ S{lz + a)] °* ~ S iM Ha "V = ^> (52) 

it is worth noting that Eq. (1521) corresponds to the Pauli equation of a 2D harmonic oscillator 

with an additional spin-orbit term. Then, it reads in coordinates representation 

(18 8 18 res 1 

<^ -—r— + -( — + ia) 2 -M 2 u 2 r 2 + 2Me + 2Muj[l-s(L z +a)y-—5(r-R) ( T 3 U = 0. (53) 

[r or or r 2 89 R J 

By using similar steps as in the relativistic case, the shifted energy levels e + u for the regular 
solution of Eq. ( |53i) are 

e + uj = u [2n + 1 + \m + a\ + s (to + a)} , with n = 0,1,2,... (54) 

If we drop out the energy related to the spin-orbit coupling term s (m + a) from Eq. ( 1341 , 
the obtained energy coincide with the result given in the first equation of the set (54) of Ref. 

For spinless particles and in absence of the AB effect (a = 0), we find the energy levels of 
the circular oscillator [TJ] shifted by constant value u. 

For the two cases when the irregular solution occurs, we find 

e + uj = to [2n + 1 — \m + a\ + s (to + a)] , (55) 

this result can be written explicitly in term of the parameter £ as 

e + uj =w[2n + l-2£], for m = —N , N > and s = -1, (56) 

e + uj = uj[2n- 1 + 2^], for m = —N — 1, N + 1 < and s = +1, (57) 
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where < £ < 1. 

It is interesting to note that if we drop out the energy related to the spin-orbit coupling 
term s (m + a) from Eq. (1551) . the obtained energy coincide with the result given in the second 
equation of the set (54) of Ref. p2], i.e. 

e + co = oo[2n + 1 - \m + a\] , (58) 

which is valid only for (m = —N , N > and s = — 1) or for (m = —N — 1, N + 1 < and 
s = +1). Eq. fl5*Bl can be also written explicitly in term of the parameter £ as 

e + u = w[2n+l-f], for m=-N , N>0 and s = -1, (59) 

e + u = u[2n + £], for m = —N — 1, N + 1 < and s = +1. (60) 

4 Conclusion 

In this paper, we have found the bound states and corresponding eigenf unctions of a 2D Dirac 
oscillator in the presence of the Aharonov-Bohm effect . The energy levels depends on the spin 
of particle and the AB magnetic flux parameter. It is shown that the energy spectrum can takes 
particular values due to the singular solution. These values (for fixed spin projection) depend 
only on the AB magnetic parameter £. Both energy levels of regular and irregular solutions 
may indicate more information on the role the magnetic vector potential. The nonrelativistic 
limit is also deduced. 

"The calculations presented here is another detection of an important role for the solutions 
of the wave equation which are singular at the origin. It is in fact very tempting to make the 
general observation that the rejection of singular wave functions basing on their singularity 
should be viewed with some suspicion" [T5] . 

Acknowledgments 

The authors thank Professor C. R. Hagen (University of Rochester, NY) for his remarks and 
comments. 



11 



References 



[i] 

[2] 

[3] 
[4] 

[5 

[6 

[7 
[8 
[9 
[10 

[11 
[12 

[13 



M. Moshinsky and A.Szczepaniak, J. Phys. A 22, (1989) L817. 

M. Moshinsky and Y. F. Smirnov, The Hamonic oscillator in Modern Physics, Harwood 
Academic Publishers, Amsterdam, 1996. 



P. Rozmej and R. Arvieu, J. Phys. A 32, (1999) 5367, arXiv: quant-ph/9903073 



M. H. Pacheco, R. R. Landim and C. A. S. Almeida, Phys. Lett. A 311 (2003) 93, arXiv: 



hep-th/0303142 



V. M. Villalba and A. R. Maggiolo , Eur. Phys. J. B 22 (2001) 31, arXiv: 
|cond-mat/0107529[ 



V. M. Villalba and R. Pino, Mod. Phys. Lett. B 17 (2003) 1331 , arXiv: |cond-mat/040 2348 
and references therein. 

Y. Aharonov and D. Bohm, Phys. Rev. 115 (1959) 485. 

Ya-Sha Yi, Tie-Zheng Qian and Zhao-Bin Su, Phys. Rev. B 55 (1997) 10632. 
C. R. Hagen, Phys. Rev. Lett. 64 (1990) 503. 

V. M. Villalba, Phys. Rev. A 49 (1994) 586, arXiv: |hep-th/9310010l 

C. R. Hagen, Int. J. Mod. Phys. A 6 (1991) 3119. 

Z. X. Wang and D. R. Duo, Special Functions, World Scientific Pub., Singapore, 1989, 
p.327. 

D. K. Park, J. Math. Phys. 36 (1995) 5453. 



[14] S. FLiigge, Practical Quantum Mechanics, 2nd ed., Springer, Berlin, 1994, pp. 107-110 (Vol. 
I)- 

[15] C. R. Hagen, Spin and the Aharonov- Bohm Effect, Preprint UR-1145 (1990). 



12 



